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Abstract 

Classically, it is well known that a single weight on a real interval leads to or- 
thogonal polynomials. In "Generalized orthogonal polynomials, discrete KP and 
Riemann-Hilbert problems", Comm. Math. Phys. 207, pp. 589-620 (1999), we 
have shown that m-periodic sequences of weights lead to "moments", polynomi- 
als defined by determinants of matrices involving these moments and 2m + 1-step 
relations between them, thus leading to 2m + 1-band matrices L. Given a Dar- 
boux transformations on L, which effect does it have on the m-periodic sequence 
of weights and on the associated polynomials ? These questions will receive a 
precise answer in this paper. The methods are based on introducing time param- 
eters in the weights, making the band matrix L evolve according to the so-called 
discrete KP hierarchy. Darboux transformations on that L translate into vertex 
operators acting on the r-function. 
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Classical situation: a weight and tridiagonal matrices. A single weight 
p(z), z e R, naturally leads to a moment matrix, 



m n = (pij) 



0<i,j<n-l 



0<M<n-l 



where (/, g) = j R fgdz and where Pj(z) := z^p(z). In turn, the moments lead to a 
sequence of monic orthogonal polynomials 





( Poo 


■ p0,n-l 


1 


\ 


1 

Pn{z) = , , det 
det m n 


Pn-1,0 


Pn—l,n—l 


z n- 


-1 




V PnO 


Pn,n—\ 


z n 


/ 



thus satisfying 



R 



p k (z)p t (z)p(z)dz = 5 ke h k . 



Then, as is classically well known, the vector p(z) = (po(z) , p\(z) , p 2 (z) , ...) of polyno- 
mials leads to tridiagonal matrices L, defined by zp(z) = Lp(z). 

Periodic sequences of weights and 2m + 1-band matrices. Instead of the 
classical situation, where Pj(z) = z^p(z), we consider an u m-periodic" sequence of 
weights p(z) := (pj(z))j>a on R; i.e., satisfying 



z m pj{z) = p j+m (z); 



in other words, 



P — (po> Pi; •••) Pm-u z m poi •••) z m p m _i, z 2m p , z 2m p m _i, 



(0.1) 



(0.2) 



This leads naturally to a 2m + 1-band matrix! Indeed, to this sequence and the inner- 
product (/, g) = J R fgdz, we associate, by analogy, the semi-infinite "moment matrix 
moo(p), where 



„(p) := (^(p))o< W <„-i : = ((^^i^)))o< i)i <n-l 



(0.3) 
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the determinant 

D n (p) := detm„(p), 
and the infinite sequence of monic polynomials, where pij = Pij(p), 



Pn{z) = . , det 
D n {p) 


( Poo 


P0,n-1 


1 


Pn-1,0 ■ 


pn—l,n—l 


z n-l 




\ PnO 


pn,n—l 


z n 



D n (p) 



det(zpi 



Pi+l,jJ0<i,j<n-l- 



(0.4) 



The second formula for p n (z) will be discussed in Lemma 2.2. Throughout the paper, 
the _D n (p)'s are assumed to be non-zero. Then the sequence p n (z) gives rise to a semi- 
infinite matrix L, defined by 

z m p(z) = Lp(z), (0.5) 
where L is a 2m + l-band matrk/Q this was established by us in M and a sketch of 



the proof will be given in Proposition 2.3. Moreover, Griinbaum and Haine |T5| had 
produced a sequence of "5-step polynomials", satisfying a fourth order differential 
equation and related to the classical Krall orthonormal polynomials. As we shall see, 
these polynomials are very special cases of our theory. We conjecture that all sequences 
of polynomials satisfying 2m + 1-step relations of the precise form (0.5) are given by 
generalized periodic sequences of weights, a slight generalization of (0.1), and limiting 
cases thereof. See definition 3.2. 

In Theorems 0.1 and 0.2, we shall be conjugating with the following matrices: 



(3A° + A 



/ A) i o o \ 

/3i 1 

o o p 2 1 

a, 

v '•• J 



and A T (3 + I 



( 1 \ 

0o 1 

ft 1 

o o p 2 1 

V '-J 



where A is the semi- infinite shift matrix A := 
the semi-infinite case, AA T — I ^ A T A. 



(8i,j-i)i,j>a, i-e., (Av) n = v n+1 . Note, in 



Theorem 0.1 (LU-Darboux transforms) The Lower-Upper Darboux transform 

L-\ m I^L- \ m I := (/3A° + A)(L - \ m I)(f3A° + A)" 1 (0.6) 
maps L into a new 2m + l-band matrix L, provided 
<WA) 



n 



$n(A) 



with arbitrary $(A) = ($ n (A)) n > E (L — A m /)" 1 (0, 0, ...). 



1 zero everywhere, except for m consecutive subdiagonals on either side of the main diagonal. 
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The null-space (L — \ m I) ^0,0, ...) is m- dimensional with basis vectors given by 

Dn(p) 



'D (p (k A 

$( fc )(A) = ( " V / ] , for 1 < k < m, 



n>0 



where 

p< k \z) := (uj k X - z)p(z) = (VA - z) Po (z), (u k X - z) Pl (z), . . .) . (0.7) 
The LU-Darboux transformation L — X m I i — > L — X m I associated with each 

Pn = 77T~~ , for fixed 1 < k < m, 

&n (X) ~ ~ 

induces a map on m-periodic weights 

p(z)^pW(z), (0.8) 

with p( k > leading to the 2m + 1-band matrix L. 

Remark : Section 5 (Theorem 5.1) contains the proof of a more general statement, 
involving linear combinations of $( fc )(A) . 

Theorem 0.2 (UL-Darboux transforms) . The Upper-Lower Darboux transform 

L - X m I ^L-X m I := (A T (3 + I)(L- X m I)(A T (3 + J)" 1 , (0.9) 
maps L into a new 2m + 1-band matrix L, provided^ 

^ = - % j 7TT mth $(A) = ($ n (A)) n>0 G(L-A m /)- 1 span(e 1 ,e 2 ,...,e m ). 

The (quasi) -null vectors $(A) of L — X m I depends projectively on 2m — 1 free parameters 
a , o m _i, 6 , &m-i Q « re fry 

$(A) = f (-lr 1 ^^! , (o.io) 



D n (p) 



n>0 



where 

p = (p , pi, pm-i, z m p , z m p m -i, z 2m p 0} z 2m p m _i, ...) , (0.11) 
2 e, := (0,...,^,0,...) e R°°. 

3 The UL-Darboux transform depends on m additional free parameters, compared to the LU trans- 
form . 



4 



Adler-van Moerbeke : Band matrix 



Nov 17, 2000 



}0, p. 5 



with 



Po{z) : 
Phiz) : 



m-l / 

(a k 5(z-uj k X) + b 
k=o V 



Pk{z) 



Pk-l[Z), 



for 1 < k < m — 1. 



mii/i 6 m _i 7^ 0, 



(0.12) 



T7ie UL-Darboux transform L — X m I i — > L — X m I induces a map on m-periodic sequence 
of weights 

P 1 — ► P, 

p leading to the 2m + 1-band matrix L. 



Corollary 0.3 An appropriate choice of and appropriate limits bk i— > oo and A t— > 
y«eW the following special Darboux transformation on the m-periodic weights 



P = (Po, Pl,---)^P = (Po, Pi, P 2 , • • .)> 



wii/i new weights 

PoO) : 
Pfc(z) : 



m— 1 



d \ „, , . p k {z) 



k=0 



dz J ' z 

for 1 < k < m — 1. 



with d m -i 7^ 0, 



Weights with (^-functions have been studied mainly by Krall and Scheffer [£1| and 
Koornwinder ||18|| , at least for the standard orthogonal polynomials. For recent expo- 
sitions on the subject, see, for instance, Andrews and Askey |PH| . Recently, they have 
been studied by Griinbaum-Haine [15 and Griinbaum-Haine-Horozov |16| . 



An integrable flow with initial m^ . We have introduced the method of inserting 
the time in the context of random matrices [|[ |], 23], where it has turned out to be 
very useful. In order to establish the results above, consider, as we did in [§, £§], the 
initial value problem, depending on two sequences of time parameters x = (xi,X2, ■■■) 
and y = (yi,y 2 , ■■■)■ 



dm c 
dx r 



A n m r 



dm c 
dy„ 



-m^A n , with initial m oo (0,0) = {{z\ Pj(z))) < ijj<oo , 

(0.13) 

where A is the customary (semi-infinite) shift matrix. As we shall establish in section 
2, imposing the condition 

K m m 00 = m 00 k Tm (0.14) 



4 the delta-function is denned in the standard way J f(z)5(X — z)dz = /(A). 
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on moment matrices m^ leads to 2m + 1-band matrices. This in turn, suggests the 
following useful reduction: given the times x, y G C°°, we define new times x,y,t & C°°, 

y = (yi,--,ym-i,o,y m+ i,...,y2m-i,o,y2m+i,-) 

t = (0, 0, t m , 0, 0, t^ra-, 0, 0, ts m , 0, ...), 

with 

tkm ■= Xkm ~ ykm for k = 1, 2, ... . (0.15) 

The point is that, letting evolve according to the variables x, y, i, will conserve the 
2m + 1-band form of L. The solution to the initial value problem (0.13) is given by the 
same moment matrix moo, as in (0.13), 

moo (p(z;x,y,t)) = ((z l ,pj(z;x,y,t))) , (0.16) 

but for weights, now depending on times x, y, i, defined as [] 

oo 

Pj (z; x, y, i) = e^r ^ r e ET £ S£ (-y) Pj+e (z), (0.17) 

£=0 

in terms of the initial condition p(z). The moments (0.16) give rise to the polynomials 
p n (z; x,y,i), as in (0.4), which, in turn, give rise to 2m + 1-band matrices L, via 
z m p = Lp. Then L satisfies the following equations^ in the time parameters (x, y, i), 

— = [(L^) + ,L], £- = [(L^)_,L], fori = l,2,..,m/i 
dL 



[(L%,L], z = l,2,.... (0.18) 

Vertex operators. In order to obtain formulae (0.7) and (0.12) for the weights, we 
consider two vertex operators, naturally associated with the integrable system (0.13) 
for 2m + 1-band matrices[] 

X 2 (A) := x(A _1 )e-^i tm ' x e^ 1 ^^7 e ^i yiA e - ^ 1 ~^^A. (0.19) 

5 The s/s denote the elementary Schur polynomials e^i * lZ = s n (t)z n . 
6 Note L 1 /™ and L 1 /" 1 are the right m th roots and Ze/t m th roots, so that: 



L i/ m = ( L i/my where L i/- m = A + ^ b k A k 

k<0 

V' m = ( L^ m Y where L x l m = c_i A' 1 + V Ct .A fc . 



fc>0 



r X (A) =diag (A , A, A 2 ,...). 
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The vertex operators (0.19) act on vectors of functions r(x,y,t) = (r n (x, y, f)) n >o- 111 
[||], we showed general linear combinations of them are the precise implementation of 
the Darboux transform (0.6) and (0.9) at the level of r-functions; see theorems 4.1 and 
4.2. Then in the end, we set (x,y,t) = (0,0,0), which yield formulae (0.7) and (0.12) 
for the new weights. 

It is well-known that the vertex operators generate Virasoro-like symmetries at the 
level of the r-functions, which translate into symmetries at the level of the "wave"- 
functions for band matrices. For the study of such symmetries, see Dickey (13], [L4| and 



[0]. For an extensive exposition on Darboux transforms, see the book by Matveev and 
Salle [0. 

Examples: 1. Tridiagonal matrices: A single weight leads to a moment matrix 
moo with Amoo = moo A and a tridiagonal matrix L; formulae (0.15) reduce to one set 
of times t := t = (ti,t 2 , ■ ■ ■)■ Equations (0.18) become the standard Toda lattice, with 
r-functions 

r n (t) = det m n (p(z)e^i . (0.20) 

The standard Toda lattice vertex operator, introduced by us in |2| and obtained from 
(0.19), 

X(t,A) = A~ 1 x(A 2 )eSi ^ A e " 2 ^i (0.21) 
has the surprising property that, given a Toda r- vector r(i) = (r , T\, . . .), the vector^ 

r(t)+cX(t,X)r(t) = (r n (t)+cX 2n ~ 2 e^ uy T n ^ 1 (t-2[X- 1 ])) n>o (0.22) 

is again a Toda r-vector. This precise operation can be implemented by a UL-Darboux 
transform, followed by a LU-Darboux transform and a limit. Note the UL-Darboux 
transform (resp. LU-Darboux transform) amounts, for a tridiagonal matrix, to a fac- 
torization of L — XI into an upper- times a lower-triangular matrix (resp. lower- times 
an upper-triangular matrix), and to multiplying the factors in the opposite order. The 
vertex operator above translates into adding a delta-function to the original weight. 
This establishes a dictionary between several points of view (explained in section 6) : 

L — X — L + L_ i — ► l! - X := L„L + i — ► U - y, = L'_L' + i — > L" — jj, := L' + L'_, 



p(z) i — > p(z) + c5(X - z) 



t + cXr (0.23) 



8 For aeC, define [a] = (f , %-,%-,...)£ C c 
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2. "Classical" polynomials, satisfying 2m + 1-step relations: Given moments 
fXi := (z\ po(z)), associated with a single weight po for standard orthogonal polynomials, 
satisfying for fixed integer m > 1, 



dz < oo, j > — m + 1, 



we define in section 7 new monic polynomials p^'(z), defined by a new moment matrix 
moo, which coincides with the old moment matrix = (Pi+j)i,j>o associated with the 
standard orthogonal polynomials, except for the first column. The pffl(z), defined by 

(detm n ) p n 1} (z) = 



L 


z 3 po(z) 







I 



m—1 



det 



^P-kd m—k—l 
k=0 
m—1 

/, V-l-kdm-k-l ~ c l 
k=0 

m—1 

^^2-kd m—k—l 

+ 2!c 2 

k=0 



Pi p 2 
P2 P3 
P3 pA 



m—1 

m—k—l 

+ (-l) m - 1 (m-l)! 

Cm— 1 Pm Pm+1 

fc=0 

m—1 

^ Pm-kdm-k-l pm+1 Pm+2 

k=0 



V 



m—1 

^ ] pn—kdm—k—1 
k=0 



Pn+1 Pn+2 



\ 



„m—l 



satisfy 2m + 1-step relations, i.e., 

z m p^ l \z) = Lp^\z), with a 2m + 1- band matrix L. 

It remains an interesting open question to find out whether such polynomials satisfy 
differential equations; on such matters, see section 7. 



1 Borel decomposition and the 2-Toda lattice 

In ^ we considered the following differential equations for the bi-infinite or semi- 
infinite matrix moo 



dm 



00 An 



8x n 



A n m r 



-mooA 



n = 1,2,..., 
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where the matrix A = (5jj_i)jj e z is the shift matrix; then (1.1) has the following 
solutions: 

m^y) = e Sr-" A " moo (0,0)e-Sr^ T " (1.2) 

in terms of some initial condition m^O, 0). In this general setup, the matrix is a 
general matrix, thus not necessarily generated by weights p. 
Consider the Borel decomposition = Sj" S 2 , for 

S\ G G_ = {lower-triangular invertible matrices, with l's on the diagonal} 
S2 G G + = {upper-triangular invertible matrices}, 

with corresponding Lie algebras g~,g+] then setting L\ := SiKS~ x , 



drrioo x dSx S 2 „_! _ dSi x 95*2 „_! 

Di— D 2 — Di — D 2 — — — D-l + — D 2 ^9- +9+ 

ox n ox n ox n ox n 

i—l _ c A n a— 1 



= S^mocV = SxA^r 1 = £? = (£?)_ + (£?)+ G + ^ 
the uniqueness of the decomposition + g + leads to 

dS ln _, ,„ n , dS 2 



Si 1 = (£?)-, = 



9x n <9x n 



Similarly setting £ 2 = S 2 A T S 2 1 , we find 



This leads to the 2-Toda equations for Si, S 2 and Ci, C 2 : 

= =f(/:?M 1i2 , -^ = ±(£2)^ (i.3) 
|| = [(£?)+, 4], ^ = [(££)_, 4], i = i,2,... (1.4) 

By 2-Toda theory 0] the problem is solved in terms of a sequence of tau-functions 



r n (x, y) = det m n (x, y), (1.5) 



with m n (x,y) defined below: 
bi-infinite case in G Z): 



m n (x,y) := (lHj(x,y))_ O0<iJ < n _ 1 , 

semi-infinite case (n > 0): 

m n (x, y) := {Hij{x, Z/)) < iij < n _ 1 , with r = 1. (1.6) 
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The two pairs of wave functions \I> = \& 2 ) and ^* = ^) defined by f\ 
^ 1 {z-x,y) = e^^S lX {z), x, y) = e^T (sj)' 1 x^ 1 ) 



y 2 (z;x,y) = e^?y^S 2 x(z), %(z; x,y) = ^ (sj) x^ 1 ) (1.7) 

satisfy 



and 

9 iTr _ f/'nA ,Tf \ 9 iTr* _ J ( fn\ \T,Ti* 



In ||23|| , with a slight notational modification ||, the wave functions are shown to 
have the following r-function representation: 

»,(.;*,„) = ( T -< J -' ar ''- g) eET-V) 

V T n (x,y) JneZ 

V T n {x,y) J neZ 

\ r n+i{x,y) J nez 

n^y) = ( ^y + ^K -ET^ z -A t (1 . 9) 



r n+ i(x,y) 



nez 



with the following bilinear identities satisfied for the wave and adjoint wave functions ^> 
and \J r *, for all m, n e Z (bi-infinite) and m, n > (semi-infinite) and x,y,x',y' G C°°: 

/ ^i n ( 2 ;;x, 2 /)^t m (z;x / ,i/ / )^— = f *2n(*; x, y)^ m (z; x', y') 7 — . (1.10) 
j2=oo Zttiz Jz=o Zmz 

The r-functions0 satisfy the following bilinear identities: 

r n (x- [z-\y)r m+x {x' + [ Z -\y')e^ z^ 1 dz 

- (t+ i(ar,y - [*])r ro (:c',y' + [z]) e £r(*-^->-™- (1.11) 



2 = 00 



2 = 



9 In this section, 

x(z) = diag (..., z , z , z , ...) in the bi-infinite case 
= diag (z , 2, ...) in the semi-infinite case. 

10 The first contour runs clockwise about a small neighborhood of z = oo, while the second runs 
counter-clockwise about z = 0. 
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they characterize the 2-Toda lattice r-functions. Note (1.7) and (1.9) yield 

(S 2 ) = diag(..., Tw+ ! (x 'f , ...) := h(x, y). (1.12) 



In [ 23 1 , the facts above are shown for the bi-infinite case; they can be carefully specialized 



to the semi-infinite case, upon setting r_j = for i = 1, 2, .... 

Consider the usual inner-product (, ) and an infinite sequence of weights p(z) = 
(p (z), px(z), ...). The moment matrix = ^oo(p(^)) will now depend on p(z). The 
following proposition will play an important role in this paper. 

Proposition 1.1 The solution to the equations 

^ = A n m oc , ^ = -m oc A T ", n=l,2,..., (1.13) 
dx n dy n 

with initial condition 

m^p^; 0,0)) = {(z\ Pj(z)))o<i,j<oo, 

is given by 

= pj{z; x, y))) , (1.14) 
where the weights pj(z;x,y) evolve as follows^ 

oo 

Pj(z;x,y) = e^i x * z 22s e (-y)p j+e (z), (1.15) 

in terms of the initial condition p(z; 0, 0) = (p (z), pi(z), ...). 
Proof : Indeed, one checks, from (1.15), 



oo 



= z k pj(z; x,y) |^ = -e^i XlZ ' se-k(-y)Pj+e(z) = ~Pj+k(z; x, y), 
oxk oyu e=k 

from which it follows that 

d d 
-p i:j (p(z;x,y)) = —{z\pj(z;x,y)) = (z l+k , pj(z; x,y)) = p i+kJ (p(z; x, y)), 



dx k dx k 

d d 
-Pij(p(z;x,y)) = — (z\p j (z : x, y)) = - (z\ p j+k (z;x,y)j = -p, M+k (p(z;x,y)), 



dy k ' dy k 

which is equivalent to (1.13). Here is an alternative way of checking this fact: since, 
from (1.14), 

(A k m 00 (p(z;x,y))) ij = (z t+k , pj(z; x,y)) and (m oc (p(2;; x, y))A T % = (z\ p j+k (z;x,y)), 
one checks 



11 The elementary Schur polynomials are denned in footnote 4; also = Si_fc. 
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e ET *-A» ( z \ p 3 ( Z] 0, 0)) < i)J < oo e- £T 



= £ s k(x)A k (z\ Pj (z; 0, 0)) 2 S ^) A 

fe=0 - ' J £=0 



oo 



(z t+k ,p H£ (z-0,0)) o ^ <oo s e (-y) 
/e^ x * zk z\f:s e (-y)p j+e (z-0,0)\ 

\ (-=0 ' 0<ij<oo 

(z\ Pj (z;x,y)) o< _ hj<oD . (1.16) 



2 Reductions of the 2-Toda Lattice 

Reduction from 2-Toda to 2m + 1-band matrices : 

For convenience, we define new vectors x,y,i G C°°, based on the vectors x, y G C°°: 

X = [Xi, X m -i, 0, X m +i, X2 m -1, 0, X2m+1, •••) 
y = (j/l, •••,2/m-l,0,y m+ l, •••,Z/2m-l,0,y 2 m+l, •••) 
t = (0, 0, t m , 0, 0, timi 0, 0, t^m, 0, ...), 

with 

ffcm := ^ifcm - 2/fem for A; = 1,2, ... . (2.1) 

Notice in this subsection, C\ and £ 2 are bi-infinite. In the next subsection, we shall 
specialize this to the semi-infinite case. 
Recall from section 1, 

f^oo = S^ 1 S 2 , C\ = SiAS^ 1 , C 2 = S 2 A. S^ 1 and r n = detm„. 

Proposition 2.1 Whenever r n (x,y) 7^ for all jigZ, the following three statements 
are equivalent: 

(i) A m moo = mooA Tm 

(ii) L™ = C™ ' , in which case L™ is a 2m + 1-band matrix. 
(Hi) £1, £ 2 , m oo and r n are functions ofx,y and i only. 

Also (i) or (ii) are invariant manifolds of the vector fields = A"m 00 , = 

-m^A 7 ™, n = l,2,...,. 
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Proof : Indeed, by the invertibility of Si and S2 under the proviso above, and re- 
membering the splitting m^ = Si 1 S 2 , we have that (i) holds if and only if 



C m = SlA m S -i = s^m^Sz 1 = Sim 00 K Ym S 2 1 = S 2 A Tm S^ 1 = 
Also note that (i) is equivalent to 



(2.2) 



= A^moo - m^A 



Tkm 



dx 



d d 
+ 



km 



dy, 



,fc = l,2,... 



kn , 



This is also tantamount to statement (iii), because the invariance of moo under d/dxk m + 
d/dykm implies the invariance of £1, £2 and r n . From the solution (1.2), if (i) holds 
at (x,y) = (0,0), it holds for all (x,y), and thus, by (2.2), if (ii) holds at (0,0), it also 
holds for all (x, y). ■ 

From Proposition 2.1, it follows that the Toda vector fields respect the band struc- 
ture of L :— L™ = C™, i.e., it is an invariant manifold of the flow. Therefore the Toda 
theory can be recast purely in terms of the 2m + 1-band matrix of the form 

L — °* A * 

—m<i<m 












a ~m+l{ — 1) •• 


• ao(-l) 


ai(-l) . 


1 


a-m(0) 


. a_i(0) 


ao(0) • 


.. a m _i(0) 1 












(2.3) 

with aj being diagonal matrices and a m = I. The vector fields below involve the ith 
powers L % / m = C\ and L l l m = C 2 of the right m th roots L l / m = Ci and left m th roots 
L l l m = £ 2 respectively; see also footnote 6. 

The m-reduced Toda lattice vector fields on L are as follows: 



dL 



dL 



^ = [(L*M) + ,L], ^ = [(ia_,I], fori = l,2,...,m^ 



dx 



dL 

F)f. 



dy 



1,2, 



Then L can be expressed in terms of a string of r-functions 

T~n ■ T n (x,y,t), 

which in the semi-infinite case will take on a very concrete form. 
Reduction from bi-infinite to semi-infinite 2-Toda : 



(2.4) 



(2.5) 
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In this section we focus on the Borel decomposition of section 1, but specifically for 
semi-infinite matrices = {Hij)ij>a, where it is unique. Remember the decomposition 
m oo = Si S 2 , where Si is lower-triangular, with l's on the diagonal and where S 2 is 
upper-triangular with h n = det(m n+1 )/ det(m n ) on the diagonal, by (1.12). Let h 
denote such a diagonal matrix. For any matrix m^, define S^rrioo) := Si and him^) : = 
h, as functions of the matrix m^. Following ||, we write the Borel decomposition, as 
follows 3 

moo = ST 1 ^ = (Sim^him^) («S(m^)) " . (2.6) 

It leads naturally to vectors of monic bi-orthogonal polynomials 

p^\z) = S{ moo ) X {z) = Six(z) and p^(z) = S(ml) X (z) = h{Sl)-\{z). (2.7) 

Upon introducing a formal inner-product ( , )o, where (%) % 1 zi)o = fiij, the polynomials 
p^'(z) and p( 2 '(z) enjoy the following orthogonality property, using (2.6): 



({p?,pf ) )o)..^Sim(h(Sj)-' 



h. 



(2.8) 



Letting the semi-infinite matrix evolve according to the differential equations (1.1), 
namely 



A n m r 



-m M A Tn , n = l,2, 



ox n dy n 
we have shown, in H, that the wave functions and ^ 2 have the following represen- 
tation in terms of the bi-orthogonal polynomials constructed from m^x, ?/) in (2.7): 

*i{z;x,y) = e^ zk p^(z;x,y) = e^ zk S lX (z) (2.9) 



n(z;x,y) = e-^ z ' k h- 1 p^(z' 1 ;x,y) = e-^-\S 2 1 ) T X (z- 1 ) 

with the p n y s being expressed in terms of r-functions r n of 2-Toda: 

,T n (x- [z~\y) _ f2 w x rjr^x.y + [z- 1 ]) 



pH\z;x,y) = z n ln ^ / ±»> , pW(z;x,y) = z 



and 



Tn(x,y) 

T n (x,y) = detm n (x,y) and h n 



r n (x,y) 

T n +l(x,y) 



(2.10) 

(2.11) 
(2.12) 



r n (x,y) 

In || , we have shown the following matrix representation for the bi-orthogonal polyno- 
mials, which then leads, using (2.7), to a representation of the lower-triangular matrices 
Sim^) and S(mL): 



(2.13) 





^ 1^00 


■ /^0,n-l 


1 


Pn\z;x,y)= det 

T n [X,y) 


/^n-1,0 


fJ"n—l,n— 1 


z n-l 






/^n,n— 1 


z n 
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( Poo 


•• Pn-1,0 


1 


P { a ) {z]x,y)= det 
r n (x,y) 


P0,n-1 


Pn—l,n—l 






V P0,n 


(J>n—l,n 





(2.14) 



Assume now the moments /i^ are given by weights = (p (2;), px(z), ...); then 
r n (x,y) = det (V, x, y))) 0<ij<n _ 1 = D n(p(x,y)), 
where pj(z;x,y) is given by (1.15), i.e., 

oo 

Pj{z;x,y) = e^i * i2 ^ s e (-y)p j+e (z). 



e=o 



Lemma 2.2 In the context of Proposition 1.1, the polynomials above have the following 
alternative representation in terms of the entries pij = (z\ pj(z;x,y)) of m: 



Pn ] {^,x,y) 



det({z\ (X- z)pj(z;x,y))) 



0<i,j<n-l 



det((^,p j (z;x,y))) <. J < n _ 1 
det(A^ij — Pi+i,j)o<i,j<n-i 

Tn(x,y) 

det({z i ,\p j (z;x,y)-p j+1 (z;x,y))) < id < n _ 1 

3<n-l 

detjXpij - Pij + i)o<ij< n -i 

T n (x,y) 



(2.15) 



(2.16) 



Proof : The proof follows from the representation (2.11) of p^ above, the repre- 
sentation (1.5) and (1.6) of r n , the representation (1.15) of pj and from the following 
identities: 

Xpij(x - [A _1 ],y) := A (z\pj(z;x - [A _1 ],y)^ 

= x(z\ e ET ( Xl '^) Z '±s e (-y)pM z )) 
\ e=o I 



z\ (A - z)p j (z;x,y) / 
= Xpij(x,y) - p i+1:j (x,y), 
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and 



\lMj(x,y + [\ *]) : = \(z\pj(z\x,y + [A 1] 

I oo 

= A (^eSr^2s < (-y-[A- 1 ])^;0,0)) 
\ e=o I 

I oo \ 

\ £=0 / 

= \fJ.ij(x,y) - fiij+i(x,y), 

which is based on the following identity: 

°° _ —i 







= A5>„H/)Wi 



A 



= J2( Xs n(-y) -Sn-i(-y)) 





3 n . 



Corollary 2.3 Given weights po, Pi, ■ ■ ■ ,p n -i, the following identity holds: 



det((z\ (A - z)pj(z))) <ij< n -i = det 



/ (z°,p (z)) ••• (z°,p n ^(z)) 
{ (z n ,p (z)) ••• (z n ,p n ^(z)) 



1 \ 



Proof : From Lemma 2.2, it follows that has two alternative expressions (2.13) 
and (2.15). Equating the two leads to the identity above. ■ 

Remark : Formula (2.15) and hence (2.13) just depend on the first formula (2.11) 
and r n = det(/%) <ij< n _i, with Pij(x,y) = {z\e^ XiZ% Pj(y,t)) . The y-dependence is 
unimportant. 

3 From ra-periodic weight sequences to 2m + 1-band 
matrices 

Given the m-periodic sequence of weights 

P = (Pj)j>o = (Po,Pi, ■ ■ .,p n -i,z m po, . . . ,z m p m -i,z 2m p , . . .,z 2m p m - u . . .), (3.1) 
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consider the initial value problem 

° = A n m 00 , = -m^A 7 "-, with initial m^^O) = ((z h Pj))o<i,, 



dXr 



dy n 



,j<ooi 



and the associated 2-Toda lattice equations 



dx 



dy„ 



(3.2) 



(3.3) 



In proposition 1.1, we gave the solution to the initial value problem (3.2) in general, 
whereas in theorem 3.1, we shall give the solution for m-periodic sequences of weights. 
This extra-structure will be important, when we deal with Darboux transforms. 

Theorem 3.1 Given the initial m-periodic weights (3.1), the systems of differential 
equations (3.2) has the following solutions with regard to the time parameters (x,y,i), 
introduced in (2.1): 



m 00 (p(z\x,y,i)) = ((z l , pj(z; x,y,t) 



0<i,j<co ' 



where 



Eoo _ r \ ~v 00 — 



J2se(-y)p j+ e(z). 
1=0 



(3.4) 
(3.5) 



is an m-periodic sequence of weights. Then the polynomials p^\ with pij := pij(p(z; x,y,t)) 
and T n (x,y,t) = det m n (p(z; x, y, t)), 





( Poo 


■ P0,n-1 


1 \ 


Pn\z;x,y,t) = det 
T n (x,y,t) 


Pn-1,0 


pn-\,n-\ 






\ PnO 


Pn,n—1 


z n j 



det(zpij — Pi+ij)o<ij< n -i 



r n [x 



,y,t) 



give rise to matrices L = , defined by z m p^ = Lp^\ such that L = is a 2m + 1- 
band matrix. The matrix C\ satisfies equations (3.3) and the 2m + 1-band matrix L the 
m-reduced Toda lattice equations (2.4). 



Proof: Since 



we have 



Pj+km — z m Pji ji k — 0, 1, 2, ... , 
= ( z \ z km pj - p j+km ) 

— (z t+km , pj) - (z\ Pj+km) 
Pi+km,j Pi,j+km 

= (A^moo-mooA^).., 

V / ij 
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and so m^ satisfies (i) of proposition 2.1 at (x,y) = (0,0) and hence for all (x,y). 
Therefore, by proposition 2.1, L := C™ is a 2m + 1 band matrix. 

From proposition 1.1, we know that the expression below for m^ is a solution of 
the initial value problem (3.2). The proof of (3.4) follows the lines of calculation (1.16). 
From there one computes 

m 00 (p(z\x,y)) 

= e^ XnAn m oo (p(z;0,0))e~^y- ATn 

\ rJV V0<m<oo 
oo oo 

= ]T s n (x)A n lz\ Pj (z; 0, 0)) ]T Sr(-ym, -Itom, ...)A Tmr e-£r ^ ATr 





oo 



= ( ^s n (x)^ +n ,^s r (-t/ m ,-y 2m , ...)p j+rm (z;0,0)\ e T,? v^ 1 

\ / 0<«,i<oo 

= ( e Er^ r ^,f;s r (-t /m ,- 2 / 2m ,...)^ m p J (^o,o)\ e -Er^ Tr 

\ / 0<i,j<oo 

= /^ )e Er^ r e Er=i^ 2fcm e-Er=i^ 2fcm p 7 ( 2 ;o,o)\ e -Er^ ATr 
= (^, e Er^ r e Er=i^ fem p.( Z ;o,o)\ e -Er^ ATr 

\ rj\ ■> i ' / 0<i,j«x> 

oo 

eEr^eEr =1 ^ m p.(^ 0, 0)) 2 s K-y)A T ' 



£=0 



e Er ^^ e Er=o ~^ km J2 s e (-y)p j+e (z; 0, 0) 



\ ^=0 / 0<ij<oo 

which establishes (3.4). The rest follows from (2.13) (see the last remark of section 2) 
and Lemma 2.2. ■ 

In the following we show that m-periodic sequences of weights lead to 2m + 1 band 
matrices, using a direct proof, thus without invoking the matrices C\ and £ 2 of 2-Toda 
theory, as in Theorem 3.1. Furthermore, it will be shown that the polynomials p^ are 
"orthogonal 1 in the sense (3.7). Consider here the slightly more general definition of 
m-periodic sequences (in comparison to (0.1)): 

Definition 3.2 Generalized m-periodic sequences of weights pi satisfy the following 
condition: for j — 0, 1, 2, . . . , 

z m pj E span{p , • • • , pm+j} and z m p j (z) = c^ m+j p m+j {z) + ..., with c jjm+j ^ 0. (3.6) 



Proposition 3.3 Given a sequence of weights p (z),pi(z), . . ., the monic polynomials 
Po(z),pi(z), . . . ,Pj(z), . . . of degree 0,1,2, .. defined by 

(pi(z), Pj (z)) = 0, 0<j<i-l (3.7) 
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are given by the same formula, as in 


Theorem 3.1, namely 






( Poo 


■ Po,n-l 


1 \ 


pM = j * det 
det m n 


pn-1,0 


Pn— l,n— 1 


z n-l 




V PnO 


Pn,n—1 


z n ) 



(3.8) 



with fiij = (z l ,pj(z)), m n = det(/ijj)o<jj< n -i- Moreover, if the pi are generalized m- 
periodic, then the polynomials (3.7) satisfy a 2m + 1-step relation; i.e., for p(z) = 
(Po(z),Pi(z),...) T , 

z m p(z) = Lp(z) (3.9) 
defines a 2m + 1 band matrix L, with m bands above and m below the diagonal. 

Proof : For < k < n — 1, the inner-product of p n (z), given by the right hand side 
of (3.8), with pk(z) automatically vanishes: 

(detm n ) (p n (z),p k (z)} = det ((p i0 , Pn, ■ ■ ■ , Pik, ■ ■ ■ , Pi,n-i, Pik)i=o,...,n) = 0. 

Furthermore, the orthogonality relation (3.7) determines the monic p„'s uniquely. To 
prove the second assertion, that L is a 2m + 1 band matrix, we proceed as follows: since 

m+j 

z m pj{z) = CjrPr(z), j = 0, 1, . . ., we have 

r=0 



m+j 



' ^ Pj ^ ' CjrPr 



for all i, j > 0, 



\ r=0 / 

m+j 

= (z t+m , Pj )-Y,C J r(z\Pr(z)) 
r=0 

m+j 

pm+i,j ^ ' Cj r pi r , 
r=0 



implying for all j > 0, 



^ pm,j ^ 
pm+l,j 

\ pm+n,j ) 



m+j 
r=0 



J'' 



/ P0,r \ 
Pl,r 

\ Pn,r j 



Therefore, by (3.8) the following determinant vanishes for arbitrary n > 0, as long as 
n — 1 > m + j, 



( Poo ■ ■ ■ Po 



DM 



det 



n-l 



\ pn,0 ■ ■ ■ pn,n-l 



Pmj ^ 



{z m p n (z),pj(z)), 



pm+n,j / 
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for all j such that < j < n — m — 1. This implies that 

z m p n {z) G {polynomials q(z) (q(z), pj(z)) = 0, for < j < n — m — 1}, 

= span{p n _ m (z) , p n -m+i (z) , • • • , } 

= span{p n _ m (z), p n _ m+ i(z), . . . , p n+m (z)}; 

the latter identity is valid, because z m p n (z) has degree n + m. Therefore L defined by 
(3.9) is 2m + 1-band, as claimed, ending the proof of Proposition 3.3. ■ 

Remark : A generalized m-periodic sequence of weights can be transformed in an 
m-periodic sequence of weights, via an invertible lower-triangular transformation of the 
Pi in the sequence p(z) = (Pj(z))j >0 ; the new sequence of weights thus obtained become 
m-periodic, i.e., 

z m Pj = z m p m+r (3.10) 

Such a transformation leaves the associated polynomials (3.8) unaffected, as is seen from 
column operations in the defining ratio of determinants in (3.8). These polynomials 
then lead to 2m + 1 band matrices L, which are thus unaffected by the lower-triangular 
operations of the p«. 

4 Darboux transformations on 2m + 1-band matrices 

The vertex operators Xj(A) := Xj(x, y, i; A), introduced in the introduction (see ||), 
play a central role in this work0 : 



Xi(A) := x(^)e Ll e ^ mi a * mie ^ XlX e ^ 

X 2 (A) := ^(A-^e-^i *- Am! e ^i ^^7 e ^i *« A V^i "^^A; (4.1) 

e.g., X 2 (A) acts on the vector r(x,y,t), as follows 

(X 2 (X)r(x,y,t)) n = e-Z?^ mi e^?y* xi \- n T n+1 (x,y- [A" 1 ],*- [A" 1 ]) 

where 

^-1 ^-(m-l) ^-(m+1) 

:= (z/i — 7-,--,y m -i -~^,y m +i — -, ) 

1 m — 1 m + 1 

i-iX- 1 } := (0,...,0,t m j 0,...,0,t 2m -— - ,0,...,0,...). 

m 2m 

The following two theorems were established in || , and will be applied in section 5 
to the concrete r n 's given by r n = det m n (p), with the p n 's as in (3.5). 

12 X (A)=diag(A°,A\A 2 ,...) 
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Theorem 4.1 (XU-Darboux transform,) Given the Toda lattice on semi-infinite 
2m + 1-band matrices, each vector $(A) in the m- dimensional null-space, i.e.f^, 

m—l 

_ ^(a fc X 1 (^ fc A))r 
$(A) = - := ^ G (Lit) - A m /) _1 (0, 0, ...) 

T T 

satisfies, as a function ofx,y,i, the following equations: 

L$ = A m $ 

_ = (i , /m)+4 , ^ = (£!*)_•, ^ = (4.2) 

/or i = 1, 2, ... noi multiples of m for the Xi and yi equations. Each $(A) determines an 
LU-Darboux transform, depending projectively on the m — l parameters a^; namely 

L — X m I i — > L — X m I := (PA + A)(L- \ m I)([3A° + A)" 1 

with 

= _ $n+l(A) . , s 

^ $n(A) ' ^ 

acis on r as 

m—l 

r f = r $ = £ (a fc X 1 (a; fc A)) r. (4.4) 

Defining := (0, 0, 1 , 0, ...) G R°°, as before, we have: 

% 

Theorem 4.2 (TJL-Darboux transform,) Given the Toda lattice on semi-infinite 
2m + 1-band matrices, the space (L — A m /) _1 spanjeo, ei, e m } is 2m- dimensional 
and thus depends projectively on 2m — 1 free parameters, i.e., 

m—l 

_ J; (^ X i(^^) + heZ? t "»» m X 2 (u k \j) r 

$(A) = — := ^ 

r r 

G (L(t) - \ m iy l span{e , e u e m }. 

T7ie vector $(A), as a function of x,y,i , satisfies the same equations (4-2) and deter- 
mines a UL-Darboux transform, with the same (3 as (4-3) (but depending projectively 
on 2m — 1 free parameters): 

L - \ m I ^L-\ m I : = (A" 1 /? + I)(L - A m /)(A" 1 /5 + J)" 1 ; 

it induces a map on r : 

m—l 

T ,_> f = A- X (r$) = A" 1 £ (afcXi^^A) + 6 fc eSr i - A " n X 2 (^ fc A)) r. 

fc=0 



13 



w is a primitive mth root of unity. 
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5 Proof of Theorems 0.1 and 0.2: Induced Darboux 
maps on m-periodic weights 

In order to prove Theorems 0.1 and 0.2, we apply Theorems 4.1 and 4.2 to the r- 
f unctions given by 

r n (x,y,t) = D n (p(z;x,y,t)) := D n (p Q (z; x,y,t), pi(z; x,y,t), . . .) = detm n (p(z; x, y, i)) 
with 

oo 

Pj (z;x,y,t) = e^T^ e EZJ^J2s £ (-y)p j+e (z), (5.1) 

1=0 

as in (3.5), where the initial condition p(z) = {pj{z)). >Q forms an m-periodic sequence 
of weights. We now perform Darboux transformations on L(x,y,t), which satisfies the 
m-reduced Toda lattice equations (2.4). Then, in the end, put x = y = i = 0. Theorems 
5.1 and 5.2 are the precise analogues of Theorems 4.1 and 4.2: 

Theorem 5.1 (XU-Darbouxj The Darboux transform for a semi-infinite 2m+l-band 
matrix, generated by the m-periodic sequences of weights p(z; x, y, i) above, 

L-\ m I^L- X m I = (PA + A)(L - \ m I)(f3A° + A)" 1 , (5.2) 

defines a new 2m + 1-band matrix L, provided (uj is a primitive mth root of unity) 

$ n+1 (A) _ ET=o ^i(w k X)D n (p(z;x,y,t)) 

^n(A) D n (p(z;x,y,t)) 

Case 1 For the special choice 

^(fe)/xx y^i{w k X)D n (p(z]x 1 y,t)) 
^ (A) = ^ D„ip{::x. t).l)) 

with arbitrary, but fixed 1 < k < n, the Darboux transformation maps T n (x,y,t) = 
D n (p) into a D n associated with a new m-periodic sequence of weights: 

D n (p(z;x,y,i))^ D n = D n (p(z; x,y,t))&*\\) 

= a k D n {{uj k \- z)p(z;x,y,t)) (5.4) 

Case 2 A general linear combination 

* / \ \ _ Y%=o ak^i{w k \)D n {p{z- x,y,t}) 
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leads to the map 

T n (x,y,t) = D n (p(z;x,y,t) i-> f n (x,y,t) = D n (p(z;x,y,t))®£\\) 



m—l 



a k D n ({uj k \ - z)p(z; x, y, t)) 

fc=0 

(-irdet((^,p ),(z i ,p 1 ),...,(^,p n )) 



0<i<n 

(5.6) 



where 



m—l 



p : = ^2 a k 5(z - u k \) 



k=0 



p e : = pe-!(z;x,y,t), for £ > 1, (5.7) 



and 



a k = a k e 



Remark : For the general case (case 2), (5.6) is the determinant a (n + 1) x (n + 1) 
matrix, instead of n x n. Therefore, to the best of our knowledge, this r-function is 
not generated in the usual way, as a determinant of the n x n upper-left corner of the 
moment matrix. If all but one of the a k s vanish, as in case 1, then the r-functions are 
generated in the usual way, as appears immediately from the second identity of (5.4). 
In the next statement, this problem will be absent. 

Theorem 5.2 (TJL - Darboux^) The Darboux transform for a semi-infinite 2m + 1- 
band matrix, arising from a m-periodic weights p(z; x, y, i), 

L — \ m I i — ► L — \ m I = (A T f3 + I)(L - \ m I)(A T (3 + (5.8) 

maps L into a new 2m + 1-band matrix L, provided (with D(p) := (D (p), Di(p), ...)), 

$ n+1 (A) (S)£o 1 (a fc X 1 (a;*A) + 6 fc eE^X 2 (^A))D(p)) n 

It acts on T n (x, y, i) = D n (p(z; x, y, t)) as follows 

r n := D n (p(z;x,y,t)) i-> f n = £> n _i(p(z;x, j/, t))$„_i(A) 

= (-l) n - 1 det((z\p ),(z\p 1 ),...(z\p, 

with 



n-1 



0<j<n-l 



Po- = Po{z;x,y,t) := ^ [a k 5{z - u A) +b k — - — — 

fc=0 \ Z A ) 

Pe- = pt-i(z;x,y,t), for £ > 1, (5.10) 
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where 

m— 1 

~a k = a k e^^Y e ^T^ m ) b k = -X m ' k ]T hj e^>^ xy ^ k . (5.11) 

3=0 

Ifb m -i 7^ ; then the p , pi, . . . form a generalized m-periodic sequence. 

Remark : Although the new sequence p(z; x, y, t) is generalized m-periodic in the 
sense of (3.6), it does not lead to a solution of the differential equations (3.2); in 
other words, it only satisfies (3.5) in the x and t variables, but not in the y variable. Of 
course, the matrix L remains a 2m + 1-band matrix, since it is effectively constructed 
from the new polynomials p n (z;x,y,t), defined by (3.8) with the new p's; see remark 
at the end of section 3. 

Corollary 5.3 An appropriate choice of a k , and appropriate limits b k i— > oo and A i — > 
in Theorem 5.2 yield the following Darboux transformation on the weights p(z;x;y;i): 

P = (Po, Pi, P2, ■ • •) >-»• P = (p"o, Pi, P2, • • •), 

where 

^ ( ( d \ k x( \ _l a Pk(z;x;y;t)\ 
Po = 2^ \°k \ J o(z) + d k — I , d m -! ^ 0, 

pi = p^i(z;x,y,t). (5.12) 

Before proving theorems 5.1 and 5.2 and corollary 5.3, we need a crucial Lemma : 

Lemma 5.4 The following two identities hold for the m-periodic sequences of weights 
of (5.1): 

Xi(A)£> n (p) 

= e^~ UmXim e^^ xl D n ((\- z)p) 

= e ^ l ^ m e^^\-l) n det({z\5{z - A)), (z\p ), (z\ p„_ 1 )) <i<n, 

(5.13) 

A- 1 eS^i^ Aim X 2 (A)L> n (p) 
= e ET^ x \-i) n - 



-i 



det[lz\ ^_ zm Pr \,(z\po),...,...{z\p n -2) 



0<i<n-l 

(5.14) 



with all the pj 's in the determinants above evaluated at x, y, t according to formula (5.1). 

Proof : : Here we use the first solution of (3.4) (and its calculation in the proof of 
Theorem 3.1), and in the second equality, we use the familiar formula e~^ u ^ % — 1 — u. 
So, one computes, using Xi(A), as defined in (4.1): 
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Xi(A)D„(/o(z;a;,j/,*)) 

£ =1 t£m* gZ^l x i A g Z^l im dt ime L^i i dxi 

• det {(<^-(z; 0,0,0)^^ 

det | (<v, Pj - (z, o, o, o) e - e~ f(i r e Er s ^ r e Er o< ^ ^ e - rr 

= e TZJ^ xtm e ^ iXi D n ({\ - z)p(z;x,y,t) 



0<i,j<n-l 



upon bringing A n in the nxn determinant, and using again the first expression (3.4) for 
moo. But, using (1.5) and (1.14), we compute, where in this calculation pi := Pi(x, y, t), 

D n ((X-z)p) = det(<V, (A - z)poj , . . . , (z\ (A - 2)p„-i)) <i<n-i 

= det( K (z\p ),...,{z\p n _ 1 ),X t ) o<i<n , using Corollary 2.3, 
= (-1)" det (<V, 5(z - A)) , (z\ p ), . . . , (z\ Pn-i)) 0<t<n , 

using the 5-function property, thus establishing the identity (5.13). 
For future use, we shall need the following easy identities : 



. I=1 — = e ™Ei W~ = ( — ^) 1/m , (5.15) 

1 CL 

and (in the exponential, one sums over i's, not multiples of m) 



e »=i = * e im 

(1 - a" 1 ) 1 /™ 

1 -a 
1 — <i m 

= i__^L(i_ a -)-i+i/- 
1 — a 

m— 1 

= ^ a l (l - a m )~ 1+1/m . (5.16) 
o 

Notice that, for any moment matrix moo defined by m-periodic weights, 



X .. X n \ X n 

I] 



in particular, using the periodicity of the sequence pj = Pj(z; 0, 0, 0), we have 



/ ij 
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Combining these two facts, we find 
/ /A T V „//A T 



A 



Now using X 2 (A), defined in (4.1) and using (3.4) for m^, compute: 
A-ieE^ ^X 2 (A) L> n (p(z, x, y, t)) 

= X rt gZ^ y! ™ Strm g ^1 r dy r 

det ( ( z\ Pi (z; 0, 0, 0)eSr=i *™^ m e E \ e" £ ^ 

\ / 0<«j'<oo 

= A 1_ "e^ 5iA< det | ^ l ,p,(z;0,0,0)e^ OO =i^(*)' m e S-i f ™ 2 '' m e E 



0<i,j<n-l 

x r z r \ 

/ 0<ij<oo 



0<i,j<n-l 



det(/^ , e srw" e s«^ , ^w) 

\\ (1 - (z/A)-)!A 



0<i,j<oo 



° V A / e -Er^ ATr using (5.15) and (5.16), 

[1 - (^J j J 0<ij<n-l 

det ( Z \e^^m eE , rZ r ^0 \r \ 



/ / \ m\ 1/m / / \m\- 

Mt)) Mi)) 



1/m ^ ^^m^l — 1/m / 

0<i,j<oo 

e -E^ ATr | , us i n g (5.17), 

J 0<i,j'<n-l 

= A 1 - n eE^ A< det|/^ c ^r *™» rro e E ^"= o1 A " 1 ~ r ^+'-( z ' °' °' °) \ 



A m - z m / 

J 0<i,j<n-l 



0<i,j<co 



m— 1 \m— 1- 



= Ae^ ffiA * det ^ ( z\ Xm _ yTn Pj+r(*. x, y, t) 

\ \ r=0 A I 0<i,j<n-l, 
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, (z\p n - 2 (z;x,y,t)) 



0<i<n-l 



The second from the last expression is a consequence of (3.4) and (3.5), according to 
the argument in the proof of Theorem 3.1 and the linearity of (3.5) with respect to the 
measures p = (po,Pi, ■■■), while the last line is obtained by replacing the jth column Cj 
by Cj — ACj_i, 2 < j < n, in the previous determinant and using the identity: 

m— 1 \m— 1— r„ m—l \m— 1— r „ „ 

EA Pj+r _ ^ A Pj+r-1 _ Pj+m-1 — * Pj-1 



r=0 A Z r=0 A Z 



A 

= -Pj-i 



Proof of Theorem 5. 1 From theorem 4.1 (the map (4.4)), and from (5.13) of Lemma 
5.4, it follows that 

m—l 



T n = D n (p) ^ Tn = a fc X 1 (cj fc A)D„(p) 

fc=0 

m—l 

= £ e^T =i ^ e EZ^ k ^a k D n ((u; k X - z)p) 

k=0 

m—l 

= Y.~ a k D n ((^^A - Z)p) 
m—l 

= (-1)" ]T ~ a k det (V , <J(z - cu fc A)>, p ), • • • , (z\ p„_ , 



0<i<n 



The expression on the right hand side of the third identity establishes the second identity 
(5.6), whereas the last identity establishes the third (5.6), ending the proof of Case 1. 
Setting all but one a k = 0, establishes (5.4) in Case 1. ■ 

Proof of Theorem 5.2 : According to Theorem 4.2 and Lemma 5.4, the UL-Darboux 
transform (5.8) with f3 n , given in (5.9) acts on r n {z\ x, y, t) := D n (p(z; x, y, t)) as follows: 



/m-l \ 



,k=0 

(-1)"- 1 det (V, £ d k 5(z - u k X)), (z\ p ), • • • , (z\ p n _ 2 >) 

V k=0 ) 0<i<n-l 
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(m—1 m-1 ( k\\m—r 
k=0 r=0 A Z 

) 0<i<n-l 

with a k as in (5.11) and b' k = b k e^ Uu)kx) \ 

/ I m—1 m—1 \m—r /m—1 \ 

(-if- 1 det U\]T M(* - **A) + £ — — - - £ fciur* p 



\m _ v 

k=0 r=0 A * \k=0 



J T I 1 



(-ir- 1 det((^,p ),(^p 1 ),---,(z i ,p n _ 1 



/ 0<i<n-l 

0<i<n-l ' 



using the new pi defined in (5.10). 

Finally, using the 5-function property in the second identity, and using p k = p k -i 
for k not a multiple of m, we prove the new sequence is generalized m-periodic: 



z m 



m—1 / _ ym _|_ (ym \ m \ 

Po = E ( a k z m 6(z - uj k X) + b k J - Xm — >-p k {z) 



k=0 \ 

X m ^UkS(z-uj k X) + b k -^^) +"Zb kPk (z) 

k=0 \ J k=0 



= X m p (z) + J2 h-ipk(z) 
k=l 

G span {po, • • • , p m }, with the condition that 6 m _i ^ 0, 



z m p k = z m p k _ x = pfc_i+ m = Pfc+ m , for k > 1, not a multiple of m, 

establishing Theorem 5.2. ■ 

Remark : As already pointed out in the remark following the statement of Theorem 
5.2, although the sequence p(x,y,t) is generalized m-periodic in the sense of definition 
3.2, it is not m-periodic in the sense of (0.1) and it only leads to a solution moo of (3.2) 
in the x and i variables, but not in y. However, since the matrix L is computed from 
the new polynomials p n (z;x,y,t) (defined in Theorem 3.1), by z m p = Lp and since 
establishing the form of p n only depended on the x-dependence of r through p(x,y,t), 
it is indeed defined by m-periodic weights. 

Proof of Corollary 5.3 : The proof follows at once from theorem 5.2 by letting A — > 0, 
and b k — > oo, and by picking appropriate a k . ■ 

Proof of Theorem 0.1, 0.2 and Corollary 0.3 : The proofs follow from setting (x, y,t) — 
(0,0,0) in Theorems 5.1, 5.2, and Corollary 5.3. 
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6 Example 1: Darboux transform for tridiagonal 
matrices 

In this section, we specialize to the case m — 1, which leads naturally to orthogonal 
polynomials, to three-step relations, and so to semi-infinite tridiagonal matrices L. The 
LU-Darboux transform on such matrices consists of decomposing the matrices L — XI as 
a product of lower- and upper-triangular matrices and multiplying them in the opposite 
order. The UL-Darboux goes the other way around. Unlike the case of bi-infinite 
matrices, the LU-Darboux map for the semi-infinite case is a unique operation, of course 
depending on the parameter A, whereas UL-Darboux depends on a free parameter a, 
besides A. 

What is the effect of this operation on weights? Theorems 5.1 and 5.2 show that 
LU-Darboux has the effect of multiplying the weight p(z) with A — z and UL-Darboux 
divides the weight by A — z, augmented by a delta-function (a/ X)5(z — A) involving the 
free parameter a. 

In ||, we have shown that, upon letting the tridiagonal, bi-infinite matrices flow 
according to the standard Toda lattice, the LU- or UL-Darboux transforms act on 
the eigenvectors as discrete Wronskians and on the r-functions as vertex operators 
especially taylored to the Toda lattice. Both transforms depend on one free (projective) 
parameter. The reduction to the semi-infinite case cuts out this freedom for the LU- 
transform, but not for the UL-transform. 

This vertex operators technology can be used very efficiently to get the results, after 
setting t = 0; in fact one can establish a dictionary between the three points of view: 
weights, vertex operators and Darboux transforms, as summarized in (0.23); the point of 
the dictionary is contained in the subsequent theorems and corollaries. The relationship 
rests on an elementary addition formula; namely, the sum of moment determinants D„ 
and -D n _i with regard to specific weights is again a moment determinant D n , but with 
respect to a new weight: 

D n {p) + cD n ^ ((A - zfp{z)) = D n (p(z) + c5(X - z)); 

this fact is not surprising, in view of the fact that if the r = (r n ) n >o is a vector of 
r-functions for the standard Toda lattice, then the following expressions 

r(t) + cK{t, A)r(t) 

forms a Toda r-vector as well, where X(£, A) is the standard Toda vertex operator, 
defined in (0.21), and acting on r as in (0.22). 

An arbitrary weight p(z) on R yields a 1-periodic sequence (p(z), zp(z), z 2 p(z), . . .) 
and a moment matrix m^, satisfying Am^ = m^A 1 (Hankel matrix). Also 

m n (p) = (fk+j(p))o<ij<n-i, D n{p) = det m n {p), with p k (p) = / z k p(z)dz, (6.1) 
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with Dq = 1. The orthogonality relations (3.7) lead to monic orthogonal polynomials 
in z of degree n 



" Jz) = dm iet 


( Po(p) 


■ Pn-l(p) 


1 \ 




Pn-l(p) ■ 


■ P2n-2(p) 


z n-l 


, with (puPjp) = 5ijh 




\ Pn(p) 


• At2n-l(p) 


z n ) 





In turn, the semi-infinite vector of polynomials p = {p n { z ))n>o leads to a semi-infinite 
tridiagonal matrix L, defined by 



zp = Lp, with L = 



I h 1 

V 



\ 



•• J 



(6.3) 



Theorem 6.1 (i) Given the weight p(z) and A G C, the eigenvector of L, corresponding 
to the eigenvalue X, 

($n(A)) n > = (p„(A)) n > = [ Dn{ %~^ P{Z)) ) >q ^( L - ^) _1 (0. °' °> -) ( 6 - 4 ) 

specifies a unique LU-Borel factorization 

/ 1 \ / fa 1 \ 



L-XI = L_L, = 



a 1 



V 



/ 



o A 



V 



/ 



with 



$ w+ l(A) _ _ 



The LU-Darboux transform 

L — A = i — > L — X = L+L-, 
induces the following map on weights p(z): 

p(z) i — > p{z){\ - z) 



(6.5) 



(6.6) 



(6.7) 



(ii) The two-dimensional eigenspace, corresponding to the eigenvalue A and with a dif- 
ferent boundary condition at n = , is given by 



(*»(A)) 



n>0 



•f/J n ((A-^))+A t+ i (fi 
A,(p) 



e^-A/r^O.O,...). (6.8) 



n>0 
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It specifies a a -dependent family of UL-Borel factorizations, 

/ a_i 1 \ / 1 



L — X — L'L' = 



ao 



V 



■■ J 



V 



L - X = L' + L'_ i — > L' — X — L'_L' + , 
inducing the following map on weights p(z) : 



(6.9) 



with the same j3 n and « n _i as in (6.5), but with $ n defined by (6.8). This defines 
UL-Darboux transforms 



(6.10) 



(6.11) 



Proof : These statements follow immediately from setting m — 1 in Theorems 0.1 
and 0.2. 



Corollary 6.2 Consider the map L i — ► L", defined by a UL-Darboux transform fol- 
lowed by a LU -transform: 

L — X — L + L^ i — >L'-\:= L_L + i — ► L' - n = L'_L' + i — > L" - \i := L' + L'_, 

where the parameter of the first UL-Darboux map is given by 



a : = 



C/J, 



(j, — A' 

then, upon taking the limit fi — > X, the map above induces a map of weights 

p(z) i — >■ p(z) + c5(X - z). 



Corollary 6.3 Concatenating m LU-Darboux transforms with parameter pi andn UL- 
Darboux transforms with rij parameters converging to Aj (n\ + ... + n r — n), induces a 
map of weights: 

Upon picking the pi appropriately, the fraction in front of p{z) in the formula above 
disappears. 
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These statements are established by letting the moment matrix flow according 
to (1.1), and then letting the associated tridiagonal matrix L flow according to the 
standard Toda lattice (remember (L n ) + denotes the strictly upper-triangular part of 
L n ) 

BT 

— = [(L n ) + ,L), n=l,2,.... (6.12) 

In the 3-reduction of 2- Toda, only one set of times t = i = (ti,t 2 , . . .) of (2.1) remain. 
The (x,y,t) evolution (3.5) of the weight p(z) reduces to the simple formula 

r ,,-,. ■ ' kz p(z), 

which was shown in a direct way in 0, for instance;in other terms, the Toda vector fields 
(6.12) linearize at the level of the weight pt{z). The deformations pt{z) of p(z) enable 
one to define t-dependent moments p k (p t (z)), associated moment matrices m n (p t (z)), 
and t-dependent monic orthogonal polynomials p n (z;t) of degree n, with L 2 -norms 

h n (t) := [ p 2 n (t, z)p(t, z)dz = (6.13) 

JR T n [t) 

The entries of the t-dependent L-matrix are expressed in terms of the r-functions 

D n (pt) = det m n (p t ) =: r n (t), (6.14) 

as follows, 

, 9 T k+ i T k ^iT k+1 

h = Tr-log and a k _ 1 = = . (6.15) 

dti T k ri 

Setting m = 1 in the vertex operators Xi(£, A) and X2(t, A) of (4.1) leads to 

Xi(t,A) := x(A)X(t,A) and X 2 (t, A) := x(A _1 )X(-t, A)A. (6.16) 

They are generating functions of symmetries of the standard Toda Lattice and act on r- 
vectors; see 0. The vertex operator X(t, A), defined in (0.21), is obtained from Xi(i, A) 
and X 2 (t, A), as follows 

X(t,X) := lim-4^(e^^ i X 2 (t,//))" 1 X 1 (t,A) = A- 1 x(A 2 )e2:^V 2 ^^ 1 ^ 

(6.17) 

has the surprising property (in view of the non-linearity of the problem) that, given a 
vector r = (r , Ti, . . .) of Toda r-functions, the new vector (see (0.22)) 

r + X(t,A)r (6.18) 

is a new vector of Toda r-functions. For connections with vertex operator algebras, see 
V. Kac 0. 

The following statements, Theorem 6.4 and Corollary 6.5 are completely parallel 
with Theorem 6.1 and Corollary 6.2. They provide a dictionary, between the three 
points of view: 
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Theorem 6.4 (i) The eigenvector 



= X 1 (t,A)r(t) E7tiXi D n ((X-z) Pt (z)) 

1 ' ) ' r{t) ~ ' I DM 



n>0 

-1/ 



G (L(t) - AI) -i (0,0,0,...) (6.19) 
induces a LU-Borel factorization, as in (6.5), with 

d 

and 

$„(t,A) <9ti Vr n+ i / 

i/ie LU-Darboux transform L(t) — A i— > L(i) — A witt new entries b n and a n , is given by 
(6.6) in terms of the new r -function: 

ri — > f = r$ = Xi(t,A)r(t). (6.21) 

(zi^ The eigenvectors 

1 (aX 1 (t,A)+ eE^Xa^A))^) 
A rit) 
^^D n ((X-z)p t (z))+D n+1 (fg) 

E (L — A/) _1 (l, 0, 0, ...) (6.22) 

induce a UL- factorization with a and [3 as in (6.20), but with $„(£, z) defined in (6.22); 
it defines a UL-Darboux transform L(t) — X i— ► L'(t) — X, as in (6.10), with new entries 
b' n and a' n , given by (6.15) in terms of the new r-function 

t i — >?' = A _1 Ar$ = A" 1 (a X x (*, A) + e^' lAl X 2 (t, A)) r(t). (6.23) 



Corollary 6.5 Consider the map L(t) i — > L"(t), defined by a UL-Darboux transform 
followed by a LU -transform, as in Corollary 6.2, with that same choice of a. It induces 
the map (6.18) at the level of Toda r-vectors: 

DM >— D n (p t (z) + ceEr ^5{X - z)) = (1 + cX(t, A)) L>„(p t ), (6.24) 

where X(t, A) is i/ie Toda lattice vertex operator (6.17). 
14 with asymptotics A) = e£ tiA 'A n (l + 0(A -1 )). 
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Instead of using Theorems 0.1 and 0.2 to establish those results, one can prove them 
directly, using the formulae in proposition 6.6 below. In this way, classical formulae 
have a natural r-function counterpart. 

Proposition 6.6 Given the weights Pt{z), the moments pi(p t (z)) and the T-junctions 
T n (t) := D n (p t ), we have the following expressions foi^\: 

• the monic orthogonal polynomials: 







/ ^0 


■ pn-1 


1 


t s 1 
PnW t )=n( ^ 

D n {pt) 


det 


Pn-1 ■ 


■■ P2n-2 


u n-l 






\ Pn 


■ p2n-\ 





D n ((u - z)pt(z)) 

D n (Pt(z)) 
nTnit-iu- 1 }) 



U 



Tnit) 



q n -i(u\t) :-- 



Pn(x;t) 



R" U — X 



p t (x)dx 



-D r 



Dn^ptlz)) n \u-Z / 

n r n {t+\u- 1 }) 



r n -i(t) 

The Christoffel-Darboux kernels: (for hi, see (6.13)) 



0<j<n 



The addition formula: 

D n (p t (z) +cS(u-z)) 



1 



D n +l(Pt 



det 



D n +i(pt) 

Tnit-lu-^-iv-^p) 



( 


1 


V 


. v n 


1 


/'() 


Pi 


pn 


u 


Pi 


P2 


■■ Pn+l 


y U n 


Pn 


pn+1 ■ 


■■ p2n 


z)pt 


») 







[uv 



T, 



n+l(t,p) 



D n {p t ) + {{u - zfp t {z] 

(l + cK(t,u))D n (p t ). 



15 Remember [a] := (a, a 2 /2, a 3 /3, 
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This last identity hinges on the addition formula: For a n x n moment matrix m n , 
the following identity holds: 

det (m n (p) + cxn(u) ® Xn(u)) = det m n (p) + cdetm n _i ((z - u) 2 p(z)^ , 

where 

Xn{u)®Xn{v) := («V)^ n . 

7 Example 2: "Classical" polynomials, satisfying 
2m + 1-step relations 

A very natural set of "classical 1 examples is to start from a weight for the standard 
orthogonal polynomials, thus corresponding to a tridiagonal matrix C\ = C 2 - Then we 
perform two consecutive Darboux transforms on the 2m + 1-diagonal matrix L = C™ = 
C™. This has the effect of mapping a 1-periodic sequence of weights to a generalized 
m-periodic sequence of weights, thus leading to 2m + 1-band matrices. Therefore, one 
is lead to a sequence of 2m + 1-step polynomials p^ derived from the "standard' ones; 
they satisfy 2m + 1-step relations, i.e., z m p^ ) = Lp^\ with 2m + 1-diagonal L, but not 
3-step relations. 

For a general m-periodic weight sequence, for appropriate choices of j3 and (3, and 
setting A = in (5.2) and (5.8), the compound map 

L^L = (/3A° + A)L(/3A° + A)" 1 i — > L = {A T (3 + I)L{A T (3 + (7.1) 

induces, according to theorems 0.1, 0.2 and corollary 0.3, the following compound map 
of weights (assuming d m _i ^ 0): 

pi — > p = (zp , zp 1} zp 2 , ...) i — > p = \ c k 5 {k) (z) + 4^rr ) , zp , zp u ...J . 



A particularly interesting case is to start with weights having the form pk(z) = 
z k p (z), where po(z) is subjected to the following condition: 



L 


z 3 p (z) 







Then the polynomials are orthogonal with respect to the weight po(z) and the map 
above becomes 



p = (z % po(z)) < i<00 i-> p = (p , p l7 p 2 , . . .) 

(7.2) 



Po{z) 1 — ,zp ,z p ,... 
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From the general theory, this new sequence is generalized m-periodic with minimal 
period m. One checks by hand, using z m 6W(z) = for < k < m - 1, that 

m— 1 

z m Po = E {c k z m 5 ( - k \z)+d n 

k=0 
m—1 

= 22 d m -k~i zTn Po( z ) 

k=0 

m 



-k-lZ" 1 k Po( z , 



J2dj-iPj- 



The new moments fx^ = (z x ,pAz)) become 

m— 1 

PiO 



m—1 



Pi 



thus defining monic polynomials p n (z 

P 
( 



( z \Po) = 2 Pi-kdm-k-i + ^2l{-l) k k\c k 5 ik 

k=0 k=0 

( z \pj) = ( z \ z3 Po) = Pi+j for j > 1, 



(detm n ) Pn\z) 



m—1 



det 



m—k—l 

k=0 
m—1 

22 Pl-kdm-k-1 — c l 
k=0 
m—1 

^2^2-kd m—k—l 

+ 2!c 2 

k=0 



Hi 
Pi 

P3 



m—1 

J2 Pm-k-ldm-k-1 + (-l) m_1 (m - l)!c m _i H r 
k=0 



P2 
P3 
P4 

Pm+1 



m—1 

^ t pm—kdm—k—1 
k=0 



pm+1 pm+2 



\ 



Pn+1 pn+2 



(7.3) 



y m— 1 



J 



m—1 

^ ] Pn—kdm—k—1 
k=0 

which satisfy 2m + 1-step relations 

z m p( 1 \z) = Lp^\z), with a 2m + 1- band matrix L. 
Because of the fact that very special cases of these polynomials have appeared in recent 



work [|15| on pentadiagonal matrices, obtained by taking squares of the classical tridiag- 
onal matrices for the Laguerre and Jacobi polynomials, we show how our polynomials 
can be specialized to those cases. Henceforth, for notational convenience, we replace ~ 
by ~ in the map (7.2). 
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Example: 5-step Laguerre polynomials. Darboux transforms for L = C\, and 

weight po(z) = z a e~ z I[ 0:OO )(z) for a > 0. 

Setting m = 2 in formule (7.2), we find the map 



with 



P= {Po{z),zp (z),z 2 p (z),...) ^ p = (p (z),p 1 (z),p 2 (z),. . .), 



p Q (z) = T{a)(cS(z) + dS'{z)) + (b+ -)p (z), with b + 0, 

z 



Pi(z) = z i p Q {z) = z a+i e z I [0:Oo) (z), i > 1, 

obtained from formula (7.2), by setting, for homogeneity considerations and without 
loss of generality, 

c = cT(a), Ci = d T(a), d = e, d\ = b. 

The moments (z l ,pj(z)) for the original sequence are given by the following expres- 
sions 

p i:j = (z\ pj) = (z\ z J p ) = r(a + i + j + 1) , 

with polynomials^ 

/ 



det m r 





(a + 1)! 


(a + 2)! 


.. 1 \ 


a+1)! 


(a + 2)! 


(a + 3)! 




a + 2)! 


(a + 3)! 


(a + 4)! 


.. ^ 2 


a + 3)! 


(a + 4)! 


(a + 5)! 


.. z 3 


a + n)\ 


(a + n + 1)! 


(a + n + 2)! . 


z n ) 



i=0 



the latter expected, the Laguerre polynomials orthogonal with regard to the 
weight po(z). 

The Darboux transformed moments pij = (z l ,pj(z)) are given by the following 
expressions 

p i0 = (z\p ) = eT(a + i) + bT(a + i + 1) + (5 ifi c - 5 i:1 d)T(a), 

Pij = (z\pj) = (z\z 3 p ) = T(a + i + j + 1) for j > 1, 

from which one computes the Darboux transformed monic polynomials 



16 a! := r(a + 1), (or)o = 1 and (a)j = a(a - I). ..(a - j + 1). 
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(det m r 



I (a-l)!e+a!6 + (ot-l)!c 
a!e + (a+l)!6-(a-l)!d 
(a + l)!e+(a + 2)!6 
(a + 2)!e+ (a + 3)! 6 



a + 1) 
a + 2) 
a + 3) 
4) 



a 



(a + 2)! 
(a + 3)! 
(a + 4)! 
(a + 5)! 



V (a + n- l)!e + + (a + n + 1)! (a + n + 2)! 



1 \ 



(7.4) 



2"/ 



The appendix to this paper gives the first four 5-step Laguerre polynomials. 

The classical Laguerre polynomials are evidently special cases of the Darboux trans- 
formed polynomials pffl 's: 



c=d=e=0, 6=1- 



It is interesting that, in an effort to find bispectral problems, Griinbaum and Haine 
had obtained special cases of these polynomials. Their method was to perform 
two explicit Darboux transforms on the explicit square L = C 2 of the 3-step relation 
£ for the Laguerre polynomials. They found, by computation, a new matrix L and 
polynomials p(z), which coincide with ours, by setting c = d = 0, e/b = a/r in (7.4), 
and hence r ^ 0. They show they are related to Laguerre by means of a differential 
equation. Indeed, given the differential equation for the Laguerre polynomials, 



d 2 d 
B = -z— — + (z - a- 1)— , with Bp n (z) 
oz 2 oz 



np n (z) 



and the operators 



P = B 







+ r and Q = B 







+ r + 1, 



dz dz 
they show that the pffl^s and p$'s are related by the following differential equations 

Pp n (z) = (n + r)p n (z) and Qp n (z) = (n + r + l)p n (z). 

Example: 5-step Jacobi polynomials. Darboux transform for L = C 2 and 
Jacobi weight^ po(z) = (2 — z) a z 13 1[o,2](z) , for a > — 1 and f3 > 0. Here the map is 
given by p i — ► p, with 



po(-z) = ^ (c 5(2) + d ^'(z)) + po(z) ye + -J , with e^0 
= z t p (z) = (2~z) a z^I m (z) for i > 1, 



(7.5) 



17 



It is more convenient to base the Jacobi weight on [0,2] rather then [-1,1]. 
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with 



r _ 2a+/m r( a + i)r(/3 + i) 



T(a + f3 + 2) 

As in the previous example, the adjustments of constants was made for homogeneity 
reasons. 

The moments for the original sequence are given by 



Hij = (z\pj) = 2 
and the Jacobi polynomials by 

£>(*)= 1 



{a + P + i+j + 



for j > 1, 



■x 



det 



det m n 



I (/3+a+l)! 

a[ 2 f3+ a +2 ( /3+1 )| 

(/3+o+2)! 

al2 f3+ a +3 (0+ 2 )! 

(/3+a+3)! 



a!2' 3 +"+ 2 Q3+1)! 

(/3+a+2)\ 
al2 l3+ a +3 (0+2)1 

(/3+o+3)! 
a!2' 3 +"+ 4 Q3+3)! 
(/3+a+4)! 



a!2 /3+ Q +3 (/3+2)! 

(/3+a+3)! 
a[ 2 /3+ Q +4 (0+ 3 )i 

(/3+a+4)! 
a!2 /3+ a +5 (0+ 4 )1 

(,3+a+5)! 



1 \ 



a!2 /3+ a +n+i (j3+n)\ a!2^+ a +"+ 2 (0+n+l)! a! 2^+"+"+ 3 (/3+n+2)! 
V (/3+a+n+l)! (/3+a+n+2)! (/3+a+n+3)! 

— E(-2) n " fc f (a + (3 + n + k) k {(3 + n) n ^ k z k 



Z n ) 



det ', 



The Darboux transformed moments are given by 

a\(P + i)\ 



(z\ p ) = 2^+ l+1 



+ + i + 1)! 

= 2 C 



e + c5j ) + (& - 



a + P + i + 1* 



V ' w (a + /? + i+ j + 1)! J ~ ' 



and the new polynomials p^ by: 
1 



det m n 



x det 



/ a!2^+»+ 1 /3!(e+c+ ^±g±^ ) 

(/M-a+1)! 

a! 2/5+«+ 2 09+1)! (e+ (/3+ °+ 2 )^- d) ) 

(/M-a+2)! 
Q !2^ + 3 (/3+2)!(e+ i^)) 



(/3+a+3)! 



a!2 /3+ Q +2 (/3+1)! 

(/3+a+2)! 
o !2 /3+a+3 (0+ 2 )j 

(/3+a+3)! 

0,12/3+0+4 (0+3)1 

(/3+a+4)! 



a!2 /3+ a +3 (/3+2)! 

(/3+a+3)! 
a!2 /3+ a +4 (0+3); 

(/3+a+4)! 
a!2) 3+ a +5 (0+4)1 

(/3+a+5)! 



~2 



«■ 2^+"+" +1 (/?+n)l (e+ 6 ^ggggj ) a | 2 /3+a+n+2 (/3+n+l)! a! 2^+"+"+3 (0+n+2)! 



(/3+a+n+l)! 



(/3+a+n+2)! 

39 
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Again in Grunbaum and Haine had considered special cases of these polyno- 
mials. Namely, the Jacobi polynomials satisfy a differential equation, 

Bp® =n(n + a l + P+l) P $;\ 

involving the differential operator 

B = z{z-2)t^\ + {{a + (3 + 2)z-2{(3+\))^-. 



Defining 



.9 d 
P = B -(z-2)— + r and Q = B + (z - 2 — + r + a + (3 + 1, 
oz oz 



they show the p^ and p^'s, for c = 0, d = 0, e/b = r/2/3 and hence r ^ 0, are related 
by the following differential equations: 

PpW = ( n 2 + (a + P)n + r)p^ 

QpW = (n 2 + (a + (3 + 2)n + a + (3 + r + l)pW. 

This paper shows that these polynomials have a determinantal representation in 
terms of moments, defined with respect to periodic sequences of weights. Moreover, the 
vertex operator technology enables one to consider general 2m + 1-band matrices. It 
remains an interesting open question to investigate the differential equations satisfied 
by the general 2m + 1-step Laguerre or Jacobi polynomials. 

8 Appendix 

The first few 5-step Laguerre polynomials are given by the following polynomials, which, 
for convenience of notation, we did not make monic; set a = a: 

Pi (z) = (e + c + ab) z — ae + d — a 2 b — ab 
p ( 2\z) = (2 e + d + a c + 2 c + a b) z 2 

-(4ae + Qe + a 2 c + 5ac + Qc + 2a 2 b + 4:ab) z 

+ (a + 2) (2ae — ad — 3d + a 2 b + abj 

p^\z) = (Qe + 2ad + 6d + a 2 c + 5ac + 6c + 2abjz 3 

-(18 ae + A8e + 3a 2 d + 21ad + 3Qd + 2a 3 c+18a 2 c + 

52ac + 48c + 6a 2 b+18ab)z 2 
+ (a + 3) (l8 a e + 24 e + a 3 c + 9 a 2 c + 26 a c + 24 c + 6 a 2 b + 12 a b) z 
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-(a + 2) (a + 3) (6 a e - a 2 d - 7a d - 12 d + 2 a 2 b + 2 a &) 

^(z) = (24e + 3a 2 d + 21ad + 36d + a 3 c + 9a 2 c + 26ac + 24c + 6a&) z 4 

-(96 a e + 360 e + 8 a 3 d + 96 a 2 d + 376 a d + 480 d + 3 a 4 c 

+ 42 a 3 c + 213 a 2 c + 462 a c + 360 c + 24 a 2 b + 96 a 6) z 3 
+3 (a + 4) (48 a e + 120 e + 2 a 3 d + 24 a 2 d + 94 a d + 120 d + a 4 c 

+ 14 a 3 c + 71 a 2 c + 154 a c + 120 c + 12 a 2 6 + 36 a &) z 2 

- (a + 3) (a + 4) (96 a e + 120 e + a 4 c + 14 a 3 c + 71 a 2 c + 154 a c 

+120 c + 24 a 2 6 + 48 a 6) z 

+ (a + 2) (a + 3) (a + 4) (24ae - a 3 d - 12a 2 d - 47ad - 60d + 6a 2 6 + 6ao) , 

etc... 
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